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Abstract 

An abstract approximation framework and convergence theory for the identification of 
thermoelastic systems is developed. Starting from an abstract operator formulation consisting 
of a coupled second order hyperbolic equation of elasticity and first order parabolic equation for 
heat conduction, well-posedness is established using linear semigroup theory in Hilbert space, and a 
class of parameter estimation problems is then defined involving mild solutions. The approximation 
framework is based upon generic Galerkin approximation of the mild solutions, and convergence 
of solutions of the resulting sequence of approximating finite dimensional parameter identification 
problems to a solution of the original infinite dimensional inverse problem is established using 
approximation results for operator semigroups. An example involving the basic equations of one 
dimensional linear thermoelasticity and a linear spline based scheme is discussed and numerical 
results indicating how our approach might be used in a study of damping mechanisms in flexible 
structures are presented. 
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1. Introduction 


In this paper we develop an abstract approximation framework and convergence theory for the 
identification of abstract linear thermoelastic systems. The inclusion of thermal effects in the 
dynamics of flexible structures has recently received an increased amount of attention as an area of 
research as a result of the problem of solar heating of large flexible spacecraft and the effect that it 
can have on the structure’s vibrational modes and the subsequent design of efficient and effective 
control laws. 

The approach we take here is somewhat different than the traditional one usually taken by other 
authors working in the area of thermoelasticity (see, for example, [4], [5], [9], [13]). Rather than 
starting with the mathematical formuation of the basic laws of continuum mechanics and thermo- 
dynamics and then linearizing to arrive at the so called basic equations of linear thermoelasticity, 
we consider an abstract operator formuation of these basic equations set in appropriately cho- 
sen infinite dimensional spaces. More precisely, we consider an abstract second order equation of 
elasticity coupled with a first order abstract parabolic equation to describe the dynamics of heat 
conduction. We then rewrite this system as an equivalent abstract first order system set in an 
appropriately constructed product Hilbert space. We establish well-posedness using results from 
linear semigroup theory (i.e, the Lumer-Phillips Theorem) and then proceed to define a class of 
inverse or parameter identification or estimation problems. 

Our approximation results are in the spirit of the treatment by Banks and Ito in [1], We develop 
an abstract approximation framework via generic Galerkin approximations to mild solutions of the 
abstract state equations. A sequence of finite dimensional approximating parameter identification 
problems result. Using an approximation and convergence result from linear semigroup theory (the 
Trotter-Kato Theorem; stability and resolvent convergence imply semigroup convergence) we are 
able to establish that under reasonable and readily verifiable assumptions, solutions to the finite 


dimensional identification problems, in some sense, approximate a solution to the original infinite 
dimensional inverse problem. 

We note that in a certain sense, the results that we present here subsume the framework developed 
in [1]. Indeed, Banks and Ito first treat first order abstract parabolic systems and then consider 
second order systems. The abstract thermoelestic systems of interest to us here consist of one 
equation of each type appropriately coupled. By zeroing out the coupling (which the framework 
we develop below can handle) in our treatment, we obtain most of the results for first and second 
order systems derived separately by Banks and Ito in [1], simultaneously. 

In order to demonstrate that in addition to being theoretically sound, our approximation frame- 
work is useful, numerically and computationally feasible, and performs well in practice (or at least, 
simulation) we apply our general approach to the basic equations of linear one dimensioned ther- 
moelasticity. This system consists of a one dimensional wave equation describing the longitudinal 
or axial vibrations of a slender rod, coupled with the heat equation together with appropriate 
boundary and initial conditions. We illustrate how a spline based scheme can be constructed so 
as to conform to our framework and indicate how it might be used in a study of dissipation or 
damping mechanisms in a flexible structure. Numerical results carried out on a Cray XMP/48 
supercomputer are presented and discussed. 

An outline of the remainder of the paper is as follows. In the second section we define the abstract 
thermoelastic system, establish its well-posedness, and state precisely the class of identification 
problems which are of interest to us. The approximation framework and convergence theory are 
developed in the third section, and the example and numerical studies are discussed in Section 4. 
Section 5 contains a brief summary of our findings and some concluding remarks. 
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2. Abstract Linear Thermoelastic Systems 


In this section we define an abstract linear thermoelastic system, and demonstrate the existence, 
uniqueness, and continuous dependence of solutions via linear semigroup theory. We then define 
an associated identification or inverse problem. 

Let Q be a metric space and let Q be a compact subset of Q. The set Q will be known as the 
admissible parameter set. For j = 1,2 let , | 'I j } rea ^ Hilbert spaces and let jVy, IHIjj 

be reflexive Banach spaces. We assume that for j — 1,2, Vj is densely and continuously embedded 
in Hj with |y>| . < /Zj |M|j G Vj. We let V* denote the continuous dual of Vj. Then with Hj as 
the pivot space, we have Vj «-*■ Hj = H* <-► V*, with Hj densely and continuously embedded in 

V*. We denote the usual operator norm on V* by ||-||j. ,J = In the usual manner, {•,•)_, is 

understood to denote both the inner product on Hj and the duality pairing on V* X Vj, for j = 1,2. 

For each q G Q we consider what we shall refer to as an abstract linear theormoelastic system 
given by 

u (t) + C(q)u(t)+ A,(q)u(t) + L(qye(t) = t> 0 (2.1) 

6(t) + A 2 (q)B(t)- L(q)u(t) = g(t-,q), t > 0, (2.2) 

u (0) = UQ ( q ) , it (0) = u 0 (q) ,0(0) = 0 o (q) , (2.3) 

where for each t > 0,u(2) G H\ and 0 (t) G H 2 . We assume that for each q G Q,C (q),A\ ( 9 ) G 

£(V! ,VC),A 2 (q)eC(V 2 ,V 2 *),L(q)eC(V l ,V 2 *),u 0 (q) £V u v 0 (q) G H\,6q (q) G H 2 , and 
f(-;q) G Li (0,T;/fi), g (•;?’) G L x (0,T; H 2 ), for some T > 0. The operator L(q)* G L (V 2 ,V*) is 
defined by 

{L(q)’rl>,<p) l = (L(q)<p,j>) 2 ,t!>eV 2 ,<peV x ,. (2.4) 

We shall also require the following further assumptions. 

(A) (Symmetry) For each q G Q the operator A x (q) is symmetric in the sense that 

(A x (q)v,i>) 1 = (Ai (<?)V>,^)x ,<^6 V x . 
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(B) (Uniform Coercivity) For j = 1,2 there exist constants £ R, independent of q £ Q, with 

/?j > 0 such that 

( A j (q)<p,<p)j + aj \tp\ 2 j > Pj ||y»|g ,<p £ Vj. 

(Note that without loss of generality we may assume that aj > 0 ,j = 1,2.) 

(C) (Uniform Boundedness) There exist positive constants = 1, 2, 6 , and independent of q 6 
for which 

< 7 j IMI j,<p e V}, j = 1,2 
IICfaMli. <*IMI, ,ve Vi, 

P(?Ml2. < HMli ,¥> € Vj. 

(Note that the final bound above implies that J[ Xr (q)* V’Uj, < p||0|| 2 G V2, as well.) 

(D) (Nonnegativity) For each 9 £ Q the operator C(q) is nonnegative in the sense that (C ( q ) <p,<p ) , > 

o,y> e 

(E) (Continuity) For <p £ V\ and ip £ V 2 the mappings q — ► A\ (q)<p,q — » (7 (9) <^,<7 — ► X(g)* ip,q —> 
M (9)^, a nd q —> L (9) 95 are continuous from Q C Q into Vi* or U 2 * (which ever is appropriate). 
Also, the mappings q — * u 0 (q),q — ► Uo (9), and 9 — > 0 O (9) are continuous from <5 C Q into 
Vi, Hi, and H 2 , respectively, as are the mappings 9 -<• f(t;q) and 9 -*■ g(t;q ) into Hi and H 2 , 
respectively, for almost every t £ [0,T]. 

(Note that the two “strong” continuity assumptions on the operators L(q) and L(q)* can be 
replaced by the single, somewhat stronger continuity assumption that the mapping 9 — + L(q) is 
continuous from Q C Q into £(V \,V 2 ) - i.e. 9 — * L(q) is continuous with respect to the uniform 
operator topology on C{V\,V 2 ).) 

(F) (Uniform Domination) There exist functions / 0 , g 0 £ Li( 0 ,T), independent of q £ Q, for which 
|/(U 9 )li < fo(t) and |9(i;9)| 2 < 9 o(t), a.e. t £ [ 0 ,T] and every 9 £ Q. 
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Remark A comment regarding our referring to the system (2.1) - (2.3) as an abstract 


thermoelastic system would be in order here. If one thinks of (2.1) as an abstract elasticity equation 
and (2.2) as an abstract heat equation, then by making some rather simple assumptions concerning 
thermal stress, and applying some basic mechanical and conservation of energy principles from 
continuum and thermal mechanics, one could justify the form of the coupling between the two 
equations (i.e. the operators L(q ) and L(q)* ). Typically however, the basic equations (both linear 
and nonlinear) of thermoelasticity are formulated from first principles using nontrivial mathematical 
formulations of the basic laws of mechanics and thermodyamics (see, for example, [4], [9], or [13]). 
In light of this, we justify our reference to (2.1) - (2.3) as an abstract thermoelastic system by 
simply stating that the thermo-mechancial systems of interest to us here, and in particular the 
basic equations of linear thermoelasticity in one and higher dimensions (see [4], [5], and [16]), can, 
via appropriate identifications, be put in this abstract form. 


We demonstrate the well-posedness of the system (2.1) - (2.3) for each q £ Q by first rewriting 
it as an equivalent first order system in an appropriate product Hilbert space and then applying 
results from linear semigroup theory. 

Let X be the Banach space defined by X = V) X //j X 11% with norm | ■ |x given by — 

1 /2 

(|MI? + IV>li + M 2 ) • For each q 6 Q let X(q) denote the Hilbert space which is set equivalent 

to X and which is endowed with the inner product (•, -} q given by 

= <^1 ( 9 ) V ’ 1 >¥’ 2) 1 

+<*1 (v^i ■> ^2 )i + + (f?i >^2)2 1 

for (<Pi,i>i,r}i) £ X,i = 1,2. We denote the norm on X(q) induced by the inner product {•,■), by 
|-| 9 and note that it is immediately clear that assumptions (B) and (C) imply that the norms [ • |x 
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and | • | 9 are equivalent, uniformly in g for g G Q. That is, there exist positive constants m and M, 
independent of q 6 Q, for which 


m 


I? - I 'lx 


<M |-L 


(2.5) 


For each q £ Q define the operator A ( q ) : Dom ( A (g)) C X (q) —* X (g) by 


A(q) = 


0/0 

-M«) —C (g) —/(g)* 

0 /(g) -A 2 (g)J 


( 2 . 6 ) 


Dom(A(q)) = {(<p,tl>,r ] )eV l xV 1 xV 2 :A(q)(<p,ip,T 1 )eX}. (2.7) 

Theorem 2.1 For each q £ Q the operator A(q ) defined in (2.6), (2.7) above is the infinitesimal 
generator of a C 0 -semigronp, {£ (t\q) : t > 0 } of bounded linear opearors on X(q) (and therefore 
X as well ). 

Proof The result will follow from, for example, Showalter [ 12 ], Theorem IV. 4. C, if we can show 
that for some u; E R the operator A ( q ) — u?I is dissipative and the operator A — A ( q ) is surjective 
for some A > u;. Toward this end, for x = ( 9 ?,^, 7 /) £ Dom(A(q )) we have 

{A(q)x,x) q = {(ij),-Ai(q)y- C {q)i> — / (g)* g , / (g) xj) - A 2 {q)'q),(‘P,ip,r})) q 
= <^i (9)V’,¥>)i +<*1 {^, 9)1 ~ {M (g)^,^)! - (C (g) ip, ip) l 
- {L (g)* r], V'), + <£ (9) g) 2 “ <^2 (?) 9, ? /) 2 

= «i ^)i - { C (g) V>, V>)i - (2I2 (g) g, g) 2 

< y Ml + y W? - A MI2 + «> W! 

< IMI? + y + «2 |g| 2 

< {<^1 (<?)¥>, v>)i + <*i Mi} + +<*2 |g| 2 

where u — max } < ^ 1 ^ 1 , ^-,q 2 |. Recall that we assumed, without loss of generality, that Oi,q 2 > 0, 
and note that when au = a 2 — 0, we have u = 0 and that the operator A(q) is dissipative. 
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In order to eliminate some rather tedious technical details, we assume for the remainder of the 


proof that a\ = c*2 = 0 (and therefore u = 0 ). The arguments to follow remain, in substance, 
essentially unchanged in the more general case of 01,02 > 0 . Let A > 0 , let (<p, ^,77) £ X and 
consider 

(A - A(q))(u,v,8) = (<p,ip,T}); 

or equivalently 

Xu — v — (p (2*8) 

Xv + A\ (q) u + C (q) v + L ( q )* 0 = ip (2.9) 

A 6 — L(q)v + A 2 (q) 9 = 77 , (2. 10) 

Solving ( 2 . 8 ) for v and then substituting into (2.9) and ( 2 . 10 ) we obtain 

(A 2 + AC (?) + A x (?)) u + L (?)* 0 = ^ + (C (?) + A) 9 (2.11) 

-£(?)«+ (1 + A - j A 1 2 (?))0 = A - 1 (?+/,(?)¥>)• ( 2 . 12 ) 

Let H = Hi X H2 be endowed with the inner product (•,•) given by 

<*,¥) = <*,,*!)! + <*2,<h) 2 , ( 2 . 13 ) 


for $ = (^1,^2) , ^ = (^1,^2) E H , and let V = V\ X V2 be endowed with the norm ||*|| given by 


11*11 = (11*1 11! + iwii) 1/2 , 


( 2 . 14 ) 


for $ = ($x, $2) E V\ Then H is a Hilbert space, V is a reflexive Banach space and V «— ► H V * 9 

with the embeddings dense and continuous. Define the operator A\ ( q ) £ C (V, V*) by 


^ ( 9 ) 


A" + AC (?) +- A\ (?) 

~M<?) 
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L(q)' ] 

(1 + A-M 2 (?))J ’ 


( 2 . 15 ) 


Then for <!> = ($i,$2) € V\ it is easily shown that 


{A x (q)*,*)>/h N^ilii + A-Va ||^ 2||2 


(2.16) 


> /3a 11*11, 

where j3\ — min {/3 X , A” 1 /^} > 0. It therefore follows that ^a(^) is an isomorphism from V onto 
V * (see, for example, Tanabe [15], Theorem 2.2.2). Recalling (2.8), (2.11) and (2.12) and setting 


(M) = (9)' 1 O’ + (C(q) + \)<p, A- 1 {r, + L{q)<p)) (2.17) 


and 

v = Xu — (2.18) 

we obtain that A — A ( q ) is a surjection, and the theorem is proved. 

For each q E Q define x (•; q) £ C ([0,T] ; X) by 

t 

x(t;q) - S(t;q)x 0 (q) + J S (t - s;q) F (s; q) ds, (2.19) 

o 

for t E [0,7] where {S(t;q) : t > 0} is the semigroup generated by A(q) y x 0 (q) = (uq (q ) , v 0 (q ) , 
0o(<?)) G A r , and F(^;^) = (0, g (/;<?)), a.e. t £ [0,7]. We shall refer to the X- valued 

function ^(*;^) — (n(*;g),t;(';9) ,^(*;^)) given by (2.19) above as the unique mild solution to the 
abstract theormoelastic system (2.1) - (2.3). 

We now proceed to define an identification problem corresponding to (2.1) - (2.3). Let Z denote 
an observation space. For i = 1,2, .. . , v and let (*, •; z) : X xQ —> R + denote a continuous 

map from X X Q C X x Q into the nonnegative real numbers. We consider the following parameter 
estimation or inverse problem. 


(ID) Given observations {zi}j- 1 C X Z at times C X [0,7] determine parameters q E Q 

i = 1 i=l 

which minimize 


J(q) = Y / $ i (x(t i -q),q-,z i ) (2.20) 

i=i 
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where for each q E Q and t{ E [0,T], x (2 : ; q) is given by (2.19). 

We do not address the question of the existence of a solution q to problem (ID) here, but rather 
defer this discussion until the next section. In the section to follow, the existence of a minimizer 
for the functional J given by (2.20) is established as a consequence of our approximation theory. 


3. Approximation and Convergence 

Using a standard Galerkin technique we construct a sequence of finite dimensional approximations 
to the abstract thermoelastic system (2.1) - (2.3). A sequence of approximating identification 
problems result. Using a well known approximation theoretic result for linear semigroups (i.e. the 
Trotter-Kato theorem) we establish that a solution to the infinite dimensional parameter estimation 
problem, (ID), exists, and moreover, that it is in some sense approximated by solutions to the finite 
dimensional approximating identification problems. 

For j = 1,2 and for each n 3 = 1,2,... let Hj J be a finite dimensional subspace of H 3 with 
H n ’ C V-, for all n,-. Let P n } : Hj -* H™ 1 denote the orthogonal projection of H 3 onto 

j j ' j j j 3 j 

computed with respect to the {•,•) . inner product. We shall require the following approximation 

assumption. 

(G) (Approximation) For j = 1,2 

lim - v?|| = 0, <peVj. 

71/ — MX) 11 J 3 

Note that the dense and continuous embedding of Vj in Hj together with assumption (G) also yield 
that lim I P?’ <p — ^1=0, <p 6 Hj, j = 1,2. 

71/— KX> 1 J 1 3 

For each q G Q we define the operators A*’ ( q ) £ £ {Hj’) J = 1.2, C ni (q) € £ (Hj ’) , and 
L n (q) 6 £ (tf" 1 , #" 2 ), n = (nx,n 2 ), using standard Galerkin approximation. More precisely, for 
j = 1,2, and <p n > £ H™‘ , we set A n -’ ( q)<p n ’ = ip n ’ G H where ip n ’ > s defined via the Riesz 
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Representation Theorem on the Hilbert space //"' as (A 3 (9) ip n > , \ n ’ )j = (i> n ’,X n ’)j, X n ’ € Hj’ ■ 
Similarly we set C n '(^)^ ni = f 1 where ^ G H j* 1 satisfies <C (9)p ni ,X ni )i = i^ 1 ,X n ')i, 
X ni G //"', and L n {q)g> n 1 = 0 " 2 where V’" 2 G //" 2 satisfies (L (q) ip ni ,x n *) 2 = {^ n2 >X” 2 )2> 
X" 2 G /7^* 2 . We define the operator T n (9)* G £ (//^ 2 , H\') to be the Hilbert space adjoint of the 
operator L n ( q ). 

We set Uq 1 (9) = P” 1 u 0 (9) , Vq 1 1 (9) = P"’^^) and 0£ 2 (9) = P? 2 #o(q), and set f ni (t;q) = 
Pj ni /(t;g) and g n2 (t;q) = P” 2 </(/;9) for almost every t e [0,T]. We then consider the finite 
dimensional system of ordinary differential equations in H n = H "* X H ^ 2 jpven by 


u n (t) + C ni (9) u n (t) + A? (q) u n (t) + L n (9)* 9 n (t) 


(3.1) 


= / ni (t;q), t> 0, 

9 n (t) + A r (9) (0 - L n (9) u n (0 = 5 n2 (t; 9) , < > 0, (3.2) 

u n (0) = u 0 n ’ (9) , u" (0) = ( 9 ) , 0 n (0) = C (9) . (3.3) 

We next proceed to rewrite (3.1) - (3.3) as an equivalent first order system. For each n'i,ri2 = 
1,2,... and n = (711,^2) l et X n = 7/” 1 x 7/” 1 X 7/" 2 be considered as a subspace of the Banach 
space X , and for each q £ Q let X" (9) = X n be considered as a subspace of the Hilbert space 
X(q). For each q G Q let A n (9) G C (A' n ) be given by 


A" ( q ) = 


0/0 

■A^(q) ~C n '(q) -L n {q)* 

0 L n (q) -A" 2 (9) 


let F n (t;q) = (0,/ ni ( t;q),g n 2 («;<?)), a.e. t G [0,T], and let x” (9) = (n” 1 (9),v 0 "‘ (?),0J 8 (g)). 
Setting x n (/) = ( u n ( t ) , u n ( t ) , 9 n (<)), we rewrite (3.1) - (3.3) as 


x n (t) = A n {q)x n (t) + F n (t;q), t> 0 


(0) = Vo (9) • 

The solution to the initial value problemn (3.4), (3.5) is given by 

t 

x n (1; 9) - 5" (t ; 9) xj (q) + J S n (t-s- 9) P" (s; 9) ds 

0 


(3.4) 

(3.5) 


(3.6) 
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fori G [0, T], where the semigroup of bounded linear operators on X n (or X n (q)), {£" ( t;q ) : t > 0}, 
is given by 


S n (<; q) = exp (tA n (</)) , t > 0 . 

We define a sequence of finite dimensional approximating identification problems as follows. 


(ID”) Given observations {z,}^ C X Z at times {tiY i=l C xJO.T], determine parameters q n G Q 
which minimize 

J n (q) = (*•;?) >9; z<) 

t=l 

where for each q € Q and U G [0, T ] , x n ( U\q ) is given by (3.6). 


The type or kind of convergence results that we are about to summarize and discuss here have 
been presented in detail in a number of places in the literature in the context of a variety of 
identification problems for various types of distributed parameter systems (see, for example [1] 
and [2]). In the context of the problems (ID) and (ID”) defined above, our general convergence 
framework and theory takes the following form. 

Using standard continuous dependence results for ordinary differential equations (see, for exam- 
ple, [8], Theorem 1.3.2) it is easily argued that the map q — ► x n (t-,q) is continuous from Q C Q 
into X n for each t 6 [0,T], and each n G Z + X Z + . It follows therefore, that the map q -»• J n (q) is 
continuous for each n. Consequently, since Q has been assumed to be a compact subset of Q, the 
existence of a solution q n to problem (ID”) for each n is assured. 

Form a directed set from Z + x Z + in the canonical way. That is, for m = (mi,m 2 ), n - 
(nj,n 2 ) £ X Z+, we say m < n if and only if mj < nj and m 2 < n 2 . Consider next the following 
proposition. 


11 


q , then 


Proposition 3.1 If { q n } is a net from the directed set Z + X Z + into Q C Q with q n — ► 
Z n (*;g n ) -> z(/;?) for each t E [0,T] where x n and z are given by (3.6) and (2.19), respectively. 


It is straight forward to argue (see [2], Chapter III, Section 1) that if Proposition 3.1 can be 
verified, then the fact that Q is a compact subset of Q implies that the net {<? n } formed from the 
solutions to the problems (ID n ), admits a convergent subnet, { q m }, whose limit, {<?}, is a solution 
to Problem (ID). Thus if Proposition 3.1 can be verified, we have established the existence of a 
solution to problem (ID), and the fact that it is in some sense approximated by the solutions to 
the finite dimensional identification problems (ID n ). 

Since Galerkin approximation was used to construct the operators A™’ (q) ,j = 1,2, C ni (<?), and 
L n (#), the same arguments used to establish that the operator A ( q ) —ul is dissipative in the proof 
of Theorem 2.1, yield that the operators A n ( q ) - a?/ are dissipative, for all n E Z + X Z + and q E Q 
with u> independent of n and q . It follows, therefore, that the semigroups {5* n (£;<?) : t > 0} are 
uniformly exponentially bounded, uniformly in n and q. That is, 

\S n (t;q)\ g <e"*, t> 0, 

or, recalling (2.5), that 

M 

\S n (t ;q )\ x t> 0. 

Thus, using assumptions (E) - (G), Proposition 3.1 will follow immediately if we can show that 

lim \S n (t;q n ) P n x - 5’(t;g)x| A - = 0 (3.7) 

for each x 6 X and t £ [0,T] whenever q n — + q, where the projection like mapping P n : X —* X n 
is given by P n = (P? 1 <p, P? 1 ip, P? 2 r}), for ( <f,ip,r 7) £ X. 
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We establish (3.7) using a version of the well-known Trotter-Kato semigroup approximation 
theorem given in [2] (Theorem II.l. 14 on page 40). According to this theorem, the convergence in 
(3.7) will follow once we have shown that for some A > u 


lim 

n— +oo 



a h ( or 1 pnx - pn ( A - A (?)) -1 ^ 


9 " 


= 0 


(3.8) 


for each iGl wherever q n — ► q. 

In light of (2.5) and assumption (G), the convergence stated in (3.8) will follow if we can argue 
that 


lim 

71— >00 


(A — A n {q n ))~ l P n x — (A - A(<7)) -1 x 



(3.9) 


for each x £ X and some A > u, wherever q n — ► q . For simplicity, we again assume that a\ = a- 2 = 0 
and therefore that u> = 0. The proof in the more general case when a\ = a 2 > 0, is essentially the 
same. 

Suppose that 

A(q))(u,v,0) = (3.10) 


and 


(A - A n (q n )) (u ni , v" , 0 n n = P" (<p, V’, 1) 


(3.11) 


= (P?'<p,P 1 n '4>, p ? 2 rj). 

We shall show that (u n ' , v n> , 8 n7 ) — » ( u,v,6 ) as n — * oo. Recall the proof of Theorem 2.1 and let 
H n denote the finite dimensional subspace of H given by H n = H x H^ 2 ■ Note that H n C V for 
all n — (ni,ri2) G Z + x Z + . For A > 0,n = (ni,n2) G Z + X Z + , and q G Q, define A™ (9) G C(H n ) 
by 


A n x (q) 


A 2 + AC” 1 (<?) + Ay 1 (q) L n (q)* 

~L n (q) (1 + A ~ l A?(q)) 


It follows that 




(3.12) 
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for $" = ($x\$£ 2 ),$ = ($" 1 ,$” 2 ) G H n , where the inner product {•,•) on H is given by (2.13) 
and the operator A x (q) G C{V,V*) is given by (2.15). From (2.16) and (3.12) we obtain that 

W (?) *”, $") > Px ||$1 2 , G H n (3.13) 


where (5\ = min(/3i,A /? 2 ) and ||-|| is the norm given in (2.14). 

From (3.10) and (3.11), as was the case in the proof of Theorem 2.1 (see (2.11) and (2.12) ), we 


find that 


and 


A\(q)(u,0) = $ = ($i,$ 2 ), 
v — Xu — 


(3.14) 


(3.15) 


v " 1 = A u ni - 

where A x (q) is given by (2.15), $ n = ($"',$£ 2 ) e H n is given by 

^ = P?'4> + (C n ' (g") + A)P"V, 

= *~ 1 (P? 2 V+L n (<i n )P? 1 <P), 

and $ = ($i,$ 2 ) € V* is given by 

= r!> + (C(q) + \)<p, 

$ 2 = A -1 (t}+ L(q)<p). 

Using assumptions (C), (E), and (G), it is not difficult to argue that the following assertions 
hold. 

0 ) 

fim ||$" — $||, = 0, 

where || -U^ denotes the usual operator norm on V*. 

(ii) The map q -*• A\ (q) $ is continuous from Q C Q into V * for each $ G V, and in particular 


n lim o ||/lA(9")®' - >1 a(?)$||, = 0, 
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whenever q n — *• q. 


(iii) There exists a positive constant p\, independent of q G Q, for which 

||Aa(?)$||, <paII®II. $e y - 

Furthermore, if we let P n : H -*• H n denote the orthogonal projection of H onto H n with respect 
to the inner product (•, •) given by (2.13), then P"$ = (P" 1 4 , x,P 2 2 $2) > $ = ($i> $2) £ H • H 
follows from assumption (G) that 

lim IIP 71 * - I'll = 0, ?eV, (3-16) 

n — *oo 

and 

lim |P"$ - $| = 0, $ G P. 

71 — KX) 

Setting $ = (u,0) and 5'” = (u ni ,0” 2 ), from (3.13) and (iii) above, we obtain 
/?*||« n - P n »|| 2 < (A" (<?")($" - P n $),$ n - P n $) 

= (A” (g n ) $ n - A x (9) $,\ $ n - P"$) 

+ (Aa (9) $ - Aa (9") $> $ n - -P"$) 

+ (A a (9") $ - A a (,") P n $, $" - P n $> 

< ||$ n - $||, 11$" - P"$|| + || Aa (9) $ - A a (9") $||, 

x ||$* _ P"$|| + || A a ( 9 n )($ - J P n $)||. ||$ n - P"$|| 

< 2 ^ {||$-$||! + „Aa(9)$-Aa(9")$||! + , 2 aII$-^$I| 2 } 

+ y ||$ n - P"$|| 2 . 

It follows therefore, from (i), (iii) and (3.16) above, that 

Urn 11$"- P"$|| 2 < 4- lim {||$ n -$||! + ||AA(9)$-AA(9 7l )$|| 2 , 
n— >0 ° /?a n- " 50 (3A7) 

+ p\ ||$ n - p n $|| 2 } = 0 . 
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The estimate (3.17) together with (3.14), (3.15), (3.16) and assumption (G) yield that 

lim (u ni , v ni , 6 712 ) — (u, v,9) 

n — ► oo 

in X and the convergence asserted in (3.9) has been established. 

Remark We note that if the admissible parameter set Q is infinite dimensional (i.e, functional), 
its discretization could also be included within our general approximation framework and conver- 
gence theory, as well. For a detailed description of how this can be accomplished, see for example, 
[3]. 


4. Examples and Numerical Results 


We consider the basic equations of one dimensional linear thermoelasticity (see [5] [16]) 


d 2 u d 3 u 

p W~ av ( + 2 ^ Ihfdi 2 


(A + 2/z) 


8 2 u 

drj 2 


89 

+ a T (3A 4- 2/i) = / 0 , 


0 < r] < £, 


t > 0 


(4.1) 


89 

8t 


8 2 9 

; <v 


d^i 

P c — - K— + 6ql t (3A + 2/i) = go. 


dtfdt 


0 < rj < £, t > 0. 


(4.2) 


With the introduction into the first equation above of the Voigt- Kelvin viscoelastic damping term, 
equations (4.1) and (4.2) describe the longitudinal, or axial, vibrations of a thin visco- thermoelastic 
rod of length £ . In the above equations, u denotes the axial displacement, 9 the absolute tempera- 
ture, p is the mass density of the rod, A and fi are the Lame parameters, c is the rod’s specific heat 
and k is its thermal conductivity. The positive constant 6 is known as the reference temperature - 
the absolute temperature of a stress free reference state for the rod. The functions / 0 and g 0 repre- 
sent, respectively, an externally applied axial force and thermal input. The nonnegative constants 
ay and a? denote the viscosity coefficient and the coefficient of thermal expansion, respectively. 
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We are interested in studying the partial differential equations (4.1), (4.2) together with initial 
conditions of the form 

rtf 

u (0, 7 ]) - u 0 (f]) , (0, rj) = n 0 ( rj ) , 6 (0, q) = 9 0 ( r }) , (4.3) 

for 0 < T) < l, and one or the other of the following two sets of boundary conditions. 

(I) Clamped ends with temperature fixed - 

u(<,0) = u(f,£) = 0, <>0 

0(i,O) = = 0. 

(II) Clamped and insulated ends - 


u(t,0) = u(t,£) = 0, f >0 


|f(t,0) = |?(M) = 0, 


t > 0. 


We note that our theory can of course handle a variety of other boundary conditions. For our 
purposes here, however, to simply illustrate the application of our general appraoch, the two given 
in (I) and (II) above will suffice. 

We consider the problem of identifying some subset of the parameters 


q = (oy,ar, A,^,k) 


(4.4) 


over some closed and bounded subset Q contained in the positive orthant of the Euclidean space 
Q = R 5 . Once again we note that our theory, or some relatively minor modification of it, can handle 
inverse problems involving the estimation of any of the parameters, input functions, or initial data 
appearing in either the differential equations (4.1), (4.2) or initial conditions (4.3). We consider 
only the subset given in (4.4) simply for the purpose of illustration. 
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In order to apply the abstract theory developed in the second and third sections above, we make 
the following definitions, identifications and assumptions. We set H\ = 1 * 2 ( 0 , £) endowed with the 
standard L2 inner product (to be denoted by*{*, *) 1 ) and set H\ = L 2 ( 0 ^£) with the weighted inner 
product 

t 

2 = J 
0 

In the case of the boundary conditions (I) we define both Vi and V2 to be the Sobolev space 
/?o(0,£) with the boundary conditions (II) we take Vi = Hq( 0 ,£) and V2 = H 1 ( 0,£). The Sobolev 
spaces Hq( 0 ,£) and //*( 0,^) are assumed to be endowed with the usual Sobolev inner products and 
induced Sobolev norms. 

For each q E Q we define the operators Aj (q) G C (Vj,V)*) ,7 = 1,2, C (q) E C (Vi , Vi*), and 
L(q)eC(V 2 ,V 2 *) by 

£ 

(A x (q)<p)W) = j^!±D<pD1>, (4.5) 

0 

£ 

(■^2 (q) <p) (VO = f —D<pDi/>, <s,il>ev 2 , (4.6) 

J PC 
0 

e 

(C{q)<p)W) = j av + VhV’GU, (4.7) 

0 

and 

€ 

(X (g) v>) (VO = - / ar g - (3A + 2 - } W, <pev lt il>ev 2 . (4.8) 

0 

Definition (2.4) then yields that the operator Z(<?)* 6 £ ( V2 , V x * ) is given by 

£ 

(L (q)* ip) (if) = J «T(3^+2/i) ^ p6 ^ 6 y 2 . (4.9) 

0 

We assume that « 0 € #o(0>^)i u o,#o € £ 2 ( 0 ,^); /o,5o € L\ (0, T; £2(0,^)) and set = 

pfo (t,v), 9 (t, V) = (t, //) for aJ m °st every G [0,7] X [0,f]. It is then a simple and straight 

forward matter to show that the assumptions (A) - (F) stated in Section 2 are satisfied. 


18 



For each N = 1,2, ... let {<p ? }^ Q denote the standard linear 5-spline (i.e. “pup-tent” or “hat”) 
functions on the interval [0,£] defined with respect to the uniform mesh, {0,i/N,2£/N, . . . 

If for n i =2,3,... we set H? 1 = span {tf 1 }^ 1 and for each n 2 = 2,3,... we set H% 2 = 
span {<p* 2 } n ^~ l in the case of boundary conditions (I) or H^ 2 = span {<p ™ 2 } in the case of 
boundary conditions (II), then Hj’ C Vj, j = 1,2, and , using standard approximation estimates 
for interpolatory splines (see [14]) together with the Schmidt inequality (see [11]) it is not difficult 
to argue that assumption (G) is satisfied (see, for example, [10], [11]). 

With the spline bases for J/* 1 and H^ 2 defined in the previous paragraph and using the definitions 
(4.5) - (4.9), it is straight forward to write the finite dimensional initial value problems (3.1) - (3.3) 
(or equivalently, (3.4), (3.5)) in matrix form. In the examples to follow we took the observation 
space Z to be Z 2 ( 0,^) and chose the functionals as 

t 

(x, q ; z) = ((«, V , 6),q;z) = J | u- z \ 2 , 

0 

i = 1,2,... ,i/. 

We used the IMSL routine ZXSSQ , an implementation of the iterative Newton’s method/Steepest 
descent hybrid Levenberg-Marquardt algorithm to solve the resulting finite dimensional nonlinear 
least squares minimization problems (ID n ). The IMSL routine IVPAG, an implementation of Gear’s 
method for stiff systems, was used in each iteration to solve the initial value problems (3.4), (3.5) 
for a given choice of q £ Q, whenever required. A composite two-point Gauss Legendre quadrature 
formula was used to numerically compute integrals whenever necessary (in inner products when 
computing orthogonal projections, for example). All codes were written in Fortran and all com- 
putations were carried out on the Cray XMP/48 at the San Diego Supercomputer Center. We 
note however, that since we were only interested in estimating constant (as opposed to functional) 
parameters, the use of the supercomputer was probably superfluous. 
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Example 4.1 In this example we seek to 1) simply illustrate that our general approach is valid 
in practice as well as in theory and that the spline based scheme that has been proposed above 
performs satisfactorily, and 2) to obtain some feel for how good of a fit can actually be obtained. 
This second goal will aid us to a certain extent in interpreting our findings in Example 4.2 to follow 
where we consider a somewhat more realistic inverse problem. 

In order to achieve the two above stated goals the following procedure was used. So called “true” 
values for the unknown parameters and functions u and 0 were chosen. Then / 0 and g 0 were 
determined so that the partial differential equations (4.1) and (4.2) were satisfied and uo,vo> and 
0 O were chosen in accordance with (4.3). In this example we considered the boundary conditions 
specified in (I). 

We chose 

p — ay — \ ^ p — aj — c — n — 6 ~ i — 1, 
u (f, rj) = (sin t) (sin i rq) , 0 < 77 < 1 , t > 0, 

and 

6 (f , 7?) = € l [if - 7}) , 0 < 77 < 1, t > 0. 

Observations were generated according to 

z{U,rj) = (sin t ,*) (sin7r rj ) , 0 < rj < 1, 

for ti = .5z, i = 1, 2, . . . ,8 = v. The discretization levels rt\ and n<i in our spline scheme were always 
chosen the same. Thus, without fear of confusion, we shall refer to their value as n — 71 1 = n 2 . 
Note that in this case, the dimension of the approximating system of ordinary differential equations 
(3.4), (3.5) is 3n - 3 

(a) In this example we seek to identify q ~ [ay^aj). As an initial guess with which to start the 
iterative optimizing search we used q° — (.25, .25). Our results are summarized in Table 4.1a 
below. 
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J n (q°) 

J n (q n ) 

CPU(m:s) 


.99045 

miWB 

3.41 

5.13xl0 -1 

0:6.34 

8 

1.0159 



2.97x10"* 

1:46.17 




3.88 

5.74X10- 1 

5:46.50 




3.94 

2.36X10" 3 

11:45.45 


Table 4.1a 

(b) In this example we consider a somewhat more challenging problem; we seek to identify the five 
parameters q = (av,o;j, As an initial guess we used q° = (.25, .25, .25, .25, .25). Our 

results are displayed in Table 4.1b below 



Oty 


A n 


K n 

J n (Q°) 

J n {q n ) 

CPU(m:s) 

4 

.89827 

1.0720 

.99102 

1.0393 

.93712 

8.25xl0 2 

2.24xl0 _1 

1:18.14 

8 

1.0090 


.94168 



9.17xlO j 

1.86x10"* 

4:56.44 




.98076 

1.0016 

.99094 

9.34x10* 

4.41xlO- J 

27:42.44 


1.0028 

1.0072 

.97233 

1.0058 

.99443 

9.41x10* 

1.32xlO- J 

40:19.79 


Table 4.1b 


Example 4.2 In this series of tests we indicate how our approach might be applied in practice, 

and how it might perform. In particular we illustrate how it may serve as an aid in the identification 

of internal dissipation (i.e. damping) mechanisms in flexible structures. 

We consider a long, slender aluminum (AL T2024 - T4) rod of length 100 in. which is clamped 

and insulated at both ends (i.e. boundary conditions (II)). We assume that it is set to vibrating 

via a mechanical input of the form 

f aoX(e,,e 3 ] (j?)sin3x<, 0 <t?<^, 0<1<T 0 

/o( ' i ’' ) = io, n<t<T, 

where xi denotes the characteristic function for the interval I C R. We assume that there is no 
externally applied thermal input - i.e. we take 

g o (t,T)) = 0, 0 < T] < C,t > 0, 
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and we assume that the system is initially at rest. That is, we take uq = vq = 0, and 0 O = 8 , 0 < 

7j < £. 

After appropriate normalization, the laboratory or experimentally measured values for the ma- 
terial parameters appearing in the system (4.1), (4.2) are given by 

£ = 1 

p = 9.82 x 1(T 2 a T = 1.29 x 10“ 3 , m 

A - 2.0635 x 10" 1 C = 5.39611 x 10" 1 1 } 

p = 1.1111 x 10" 1 k = 7.02176 x 10' 7 . 

There is no generally accepted means to explicity measure the viscosity coefficient ay in the labo- 
ratory. Typically its value is determined via an identification procedure such as the one which we 
have developed here. We took 6 = 6.8 X lO 1 . 

With the boundary conditions given in (II), the natural mode shapes for the homogeneous sys- 
tem corresponding to (4.1), (4.2) are of the form 4>o = ( u o,9o) = (0,1), $m,j = (u m ,0 m|<7 *) = 
(sin cos ^p 2 ), m — 1,2,... j = 1,2,3, where for each m — 1,2,..,, and j = 1,2,3, 

c m j is a complex constant. Thus when a Galerkin scheme using M T 1 elements of the form 
zz (sin cos ^y 1 ) ,m = 0,1,2,... ,A/, is employed to discretize (4.1), (4.2), the resulting 
finite dimensional equivalent first order system consists of M + 1 decoupled first order systems, one 
of dimension one and the remaining M each of dimension three. We used such a scheme with the 
true values of the parameters and M = 24 to generate simulated observational data upon which 
to base our fits. Once again in our spline scheme we use n — n\ — n<i to denote the level of 
discretization. 

The presence of either, and certainly both, the Voigt-Kelvin viscoelastic term (i.e. ay > 0) 
and the thermoelastic term (i.e. a? > 0) yield some form of mechanical energy dissipation. In 
fact, when the boundary conditions are as given in (II), with a simple change of variable and the 
construction of an appropriate Lyapunov functional (or, equivalent norm), it can be shown that 
either the Voigt-Kelvin damping or the thermoelastic coupling leads to a uniformly exponentially 
stable open-loop or unforced system (i.e. §(t;q) < Mexp(u>t) ,t > 0, with cj < 0, where 
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|5(f;g) : t > o| is the semigroup corresponding to the transformed system (see [6] and [7]). The 
simple change of variable consists of transforming the absolute temperature 0 to the temperature 
variation from the time invariant average temperature over the length of the rod, 0 AVE • That such 
a constant average temperature exists, independent of the input /o and initial conditions (4.3), is a 
direct consequence of the insulated ends. Indeed, integrating equation (4.2) from 0 to l , integrating 
by parts and then imposing the boundary conditions (II), yields the desired conclusion. We note 
that even without the change of variable, the mechanical vibrations of the unforced system still 
tend to zero exponentialy fast. (In this case the absolute temperature tends to 0 AVE exponential 
fast.) Thus the presence of either the viscoelastic (Voigt-Kelvin) or the thermoelastic effects will 
exponentially damp unforced mechanical vibrations. In the series of examples to follow we illustrate 
how our scheme and general approach might be used to aid in the identification of the mechanism 
for, and/or an appropriate model for an observed dissipation of mechanical energy in a flexible 
structure. 


(a) In this test we took a 0 = Mi = .4, 5 2 = .435, T 0 = .5 and T = 5. We set a v = 0 (i.e. 
no Voigt-Kelvin damping), generated observations at times t{ = .25 i,i = 0,1,2,... ,20 via the 
modal scheme described above, and used our scheme to see if we could estimate Oct • All other 
parameters were taken as given in (4.10) and ay was held fixed at zero. For an initial guess 
we set a^ = 10“ 2 . Our results are summarized in Table 4.2a below. It is clear from the table 
that n had to be taken quite large in order to obtain an accurate estimate for ax- It is likely 
that the reason for this is that thermoelastic effects in metals are relatively slight. Consequently 
we would expect that the thermo-mechanical coupling as described by the parameter a T is not 
readily discernable by our approximation and identification procedures. 
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a? 

J n (q°) 

J n (q n ) 

' CPU(m:s) 

4 

3.0409xl0 -3 


1.187xl0 _1 

0:4.94 

8 

4.4408xl0 -3 



1:11.12 

mm 

2.3034xl0 _:J 

mmmm 


3:15.86 

mm 

BTtntnTBBI 

i 

BESlSXififl 

21:32.26 


1 

4.471x10-'" 

EUBlS&Sfli 

45:36.64 


Table 4.2a 


(b) This test is similar to the previous one described in (a) above, however, this time we introduced 
some Voigt-Kelvin damping into the system. We set ay = 10~ 2 when the observations were 
generated. This time we also took uq = 5.0 , T = 3.0, and once again set a® = 10” 2 . Our 
findings are reported in Table 4.2b. In general, the presence of the Voigt-Kelvin damping tended 
to make our task somewhat simpler and to improve the overall performance of our scheme. This 
is not surprising since Voigt-Kelvin is a rather strong form of dissipation and its presence would 
be likely to significantly increase the stability and accuracy of our numerical integrator. 


■ 

Qj 

J n (q°) 

J n {q n ) 

CPU(m:s) 

El 

2.1998xlQ- 4 

1.607xl0 _1 

1.570xl0 _1 

0:54.21 

8 

■iH.-T.ia pin— 

6.745x10"" 

6.376x10"" 

1:3.82 

16 


1.453x10-" 

9.691xl0 -a 

4:35.93 

20 



:■ 1.361x10"* 

14:17.86 


Table 4.2b 


(c) In this example we first attempt to fit a model with only thermoelastic damping to observations 
of a system which has both visco- and theoremoelastic dissipation. We then try to fit a model 
with viscoelastic damping only to a system in which the dissipation is due only to thermal effects. 
In both tests the input excitation and observations are as they were in test (b) above. 

For the first test we generated observations with ay — 10“ 2 and aj as given in (4.10). We 
then attempted to estimate the parameter ay with aj held fixed at zero. The initial guess for 
ay was taken to be ay = 2.5 x 10” 1 . Our results are given in Table 4.2c.l. 

In the second test, observations were generated with ay = 0 and aj as given in (4.10). We 
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then again attempted to estimate ay with ar fixed at zero. The initial guess was ay = 10 2 
The results for this test are given in Table 4. 2c. 2 


n 

ay 

J n (q°) 

J n (q n ) 

CPU(m:s) 

4 

7.1179xl0“ 3 

4.234xl0 -1 

1. 466x10“* 

0:10.23 

8 

1.2178x10“" 

4.122x10“* 

4.863x10“" 

0:44.62 

16 

9.7023xl0“ 3 

4.157xl0 -1 

1.064x10“! 

9:12.97 

20 

9.7812xl0“ 3 

4. 165xl0 _1 

1.135x10“" 

16:17.10 


Table 4.2c.l 


n 

ay 




4 

1.4501xl0“ 3 

3.329x10“* 

3.019x10“* 

0:8.68 

8 


msssam 

BaSaliBH 

0:24.14 

16 


KiiiSIiBfli 

kksuHI 

3:43.28 

20 


KSUSSiHi 

1 

12:53.48 


Table 4.2c.2 

For the system that we have considered here, the results given in Tables 4.2c. 1 and 4. 2c. 2 seem 
to contribute to the conclusion that the effects of thermoelastic damping are insignificant. Indeed, 
in the first test, a relatively accurate estimate for the viscosity coefficient ay can be obtained with 
a model that ignores all thermoelastic effects. The second test, although not conclusively, appeared 
to indicate that the presence of thermoelastic effects in the data will have only a minimal effect on 
the estimates obtained for the viscosity term ay 

We note at this point however, with a word of caution, that the interpretation given in the pre- 
vious paragraph is not intended as a sweeping generalization and conclusion about thermoelastic 
effects in the vibration of flexible structures. Indeed, we present these findings simply as an illus- 
tration of how our general approach and schemes could be used in such a damping study. To draw 
such conclusions would necessarily require an exhaustive investigation involving both simulation 
and computational studies as well as laboratory experimentation. 
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5. Summary and Concluding Remarks 

In this paper we have considered inverse or parameter identification problems for thermoelastic 
systems. We defined a class of abstract infinite dimensional systems consisting of a second order 
elasticity equation coupled with a first order abstract parabolic equation. We established the 
well-posedness of this system by rewriting it as an equivalent abstract first order system set in an 
appropriate product Hilbert space and then applying results from linear semigroup theory. A class of 
identification problems in which the state constraints were given by the mild solution to the abstract 
thermoelastic system were defined. We then proceeded to develop an abstract approximation 
framework based upon generic Galerkin approximation of the mild solutions. Using approximation 
results from linear semigroup theory, we were, under appropriate assumptions, able to establish that 
the solutions to the resulting finite dimensional approximating parameter estimation problems, in 
some sense, approximate a solution to the original infinte dimensional identification problem. We 
demonstrated that our approach could be applied to the basic equations of linear thermoelasticity 
and provided computational results as evidence that our schemes are not only theoretically well 
founded, but numerically sound and pactical as well. 

Many interesting and important open questions related to the identification and control of ab- 
stract thermoelastic systems remain. In particular, nonlinear systems and especially ones in which 
the material parameters are permitted to be temperature dependent need to be considered. Also, 
numerical studies using actual experimental data rather than simulation data should be carried 
out and reported on. With regard to the optimal control problem, a complete understanding of 
the asymptotic behavior of the spectrum of the operator A(q) given by (2.6) would be extremely 
useful as would it be for just the basic equations of linear themoelasticity (4.1), (4.2). Indeed, 
for the boundary equations (II) it can be shown that the eigenvalues asymptotically approach a 
straight line in the left half plane. For other standard sets of boundary conditions, in particular 
the Dirichlet/Dirichlet case given in (I), this author knows of no similar types of theoretical results. 
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Related to this is the question of uniform exponential stability of the open-loop semigroups. For 
some sets of boundary conditions, this can be established (see [7]). However, the question in general 
is still open. Finally, when doing approximation for LQG control, among the conditions required 
to establish convergence of the optimal closed-loop feedback gains, is the preservation of uniform 
exponential stability under approximation (see, for example, [6]). Even for a linear spline based 
Galerkin scheme for the basic equations of linear thermoelasticity with the boundary conditions 
given in (II), this condition does not appear to hold. Schemes for which this condition can be 
verified (other than modal schemes of course, for which it is trivially true) need to be developed 
and studied. 
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